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ABSTRACT 

o 

T-H Context. The broadband SEDs of blazars show two distinct components which in leptonic models are associated with 

C*^ synchrotron and SSC emission of highly relativistic electrons. In some sources the SSC component dominates the 

^N| synchrotron peak by one or more orders of magnitude implying that the electrons mainly cool by inverse Compton 

I collisions with their self-made synchrotron photons. Therefore, the linear synchrotron loss of electrons, which is normally 

^ invoked in emission models, has to be replaced by a nonlinear loss rate depending on an energy integral of the electron 

^~~i distribution. This modified electron cooling changes significantly the emerging radiation spectra. 

I Aims. It is the purpose of this work to apply this new cooling scenario to relativistic power-law distributed electrons, 
which are injected instantaneously into the jet. 

I I Methods. We will first solve the differential equation of the volume-averaged differential number density of the electrons, 

nj and then discuss their temporal evolution. Since any non-linear cooling will turn into linear cooling after some time, we 

HH also calculated the electron number density for a combined cooling scenario consisting of both the linear and non-linear 

HH cooling. 

(-H For all cases, we will also calculate analytically the emerging optically thin synchrotron fluence spectrum which will be 

O , compared to a numerical solution. 

I Results. For small normalized frequencies / < 1 the fluence spectra show constant spectral indices. We find for linear 

O cooling asYN = 1/2, and for non-linear cooling assc = 3/2. In the combined cooling scenario we obtain for the small 

ri injection parameter /?i — 1/2, and for the large injection parameter /32 ~ 3/2, which becomes /?i — 1/2 for very small 

C/3 frequencies, again. These identical behaviors, as compared to the existing calculations for monoenergetically injected 

Cw electrons, prove that the spectral behavior of the total synchrotron fiuencc is independent from the functional form of 
the energy injection spectrum. 
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1. Introduction tively beamed in the forward direction. Thus, if an AGN is 

(^ . . . observed straight down the jet, it is extremely bright and 

^^ Amongst the brightest sources visible throughout the ob- jg (.oHed a blazar 
r servable universe are blazars, the most extreme class of 
^^ Active Galactic Nuclei (AGN). AGN are galaxies whose 

O central super massive black hole (SMBH) is fed by tremen- xhe observed photon energy spectrum of a blazar is 
. . dous amounts of gas, which accumulate in an accretion dominated by two broad components of nonthermal radia- 
> disk surrounding the SMBH. Magnetic fields trapped in tio^ j^ leptonic models (for a review see Bottcher (2007)) 
the plasma are swirled around by the rotating disk form- the low-energetic component is attributed to synchrotron 
ing two narrow channels, commonly known as jets, along radiation of highly relativistic electrons, while the high- 
^ the rotational axis of the SMBH. Although the launching energetic component is due to inverse Compton collisions 
process in all its details is not yet fully understood (Spruit of the electrons with ambient radiation fields. Such radi- 
(2010)), it is well established that it is an important process ation fields could be of external nature (e. g. Dermer & 
to remove angular momentum from the disk. The angular Schhckeiser (1993)), like the cosmic microwave background, 
momentum is carried away by particles moving at highly radiation directly from the accretion disk, or photons from 
relativistic speeds through the jets. the accretion disk scattered in the infrared torus surround- 
In most jet models electrons (negatrons and positrons) ing the disk or in the Broad- or Narrow-Line-Regions cir- 
form the particle content of the jet, however heavier cHng the SMBH. However, since the electrons will in any 
hadronic components might also be present. The charged case produce synchrotron radiation due to their interaction 
particles are subject to several radiation processes, and due with the magnetic field of the jet, it is unavoidable that they 
to the relativistic motions the emerging radiation is effec- Compton-scatter their self-made synchrotron photons up to 

7-ray energies. This is called the synchrotron self-Compton 
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Another very important topic concerning blazars is the 
short time variaility of their emission at practically all pho- 
ton energies. Blazars have shown variabilities on the or- 
der of minutes (Aharonian et al. (2007), Ghisellini et al. 
(2009a)), and such rapid variabilities have to be explained 
by radiation and emission models. In this respect, some an- 
alytical work has been done recently regarding the possibil- 
ity of very rapid nonlinear radiation processes, which could 
be at work in leptonic radiation models of jets. Schlickeiser 
& Lerche (2007,2008) and Zacharias & Schlickeiser (2009) 
discussed a nonlinear synchrotron model which assumed 
equipartition between the magnetic field and the electron 
energy density. Schlickeiser (2009, hereafter referred to as 
RS) investigated as additional nonlinear cooling process the 
synchrotron self-Compton radiation process, which is de- 
scribed above. 

In many Blazars, Uke PKS 0048-071, PKS 0202-17, PKS 
0528+134, 3C 454.3 (Ghisellini, Tavecchio & Ghirlanda 
(2009)), S5 0014+813 (Ghisellini et al. (2009b)), and oth- 
ers from the Fermi blazar survey (Abdo et al. 2010), the 
Compton peak dominates by at least one order of mag- 
nitude over the synchrotron peak, which is even more 
pronounced when the 7-ray absorbtion of the extragalac- 
tic background light (EBL) is factored in (e.g. Venters, 
PavHdou & Reyes (2009)). 

Assuming that the Compton peak is mainly a result of SSC 
radiation one finds for the ratio of the luminosities of the 
peaks (Schlickeiser, Bottcher & Menzler (2010), hereafter 
referred to as SBM) 



^ssc 



I dv J- 



L*sYN I dV J, 



d7n(7)|7ssc| 



(1) 



Since SSC cooling relies on the synchrotron photons of the 
same electrons, the electron density n('-f) is the same in 
both cases. Similarly, the Doppler boosting factors are also 
identical (Dermer & Schlickeiser (2002)), which means that 
the ratio of the luminosities directly refiects the ratio of 
the cooling factors \ji\ (with i € {SSC, SYN}). Hence, a 
dominance of the Compton peak over the synchrotron peak 
implies that the electrons mainly cool by inverse Compton 
collisions with their self-made photons. In this case, the 
linear synchrotron cooling has to be replaced by another 
cooling process dealing with SSC radiation, which has been 
considered by RS. 

It is important to notice that the nonlinear SSC cool- 
ing only operates in the Thomson regime and does not 
deal with higher order SSC collisions, which are possible. 
Taking higher oder interactions and/or the Klein-Nishina 
(KN) regime into account, would lead to a much more com- 
plicated approach. If KN effects are important in the begin- 
ning, our approach would be at least a good approximation 
for later times, since the electrons cool significantly over 
time and will sooner or later reach the Thomson scatter- 
ing regime. However, blazar spectral energy distributions 
(SEDs) normally do not show a third broad component, 
which would be an indicator of higher order scattering or 
scattering in the KN domain. This is a hint for the valid- 
ity of the approach. We notice that this restriction to the 
Thomson regime requires according to RS that the maxi- 
mum Lorentz factor is not larger than jmax ~ 1.94-10*6^^'"^ 
(6 is the normalized magnetic field strength in units of 
Gauss) . 



RS applied the nonlinear SSC model to the simple but 
illustrative case of mono-energetic electrons that are instan- 
taneously injected into the jet. It is the purpose of this work 
to apply this new model to a more realistic scenario. We 
assume an instantaneous injection of power-law distributed 
electrons of the form Q{'f,t) oc 7^" 5 [t). This describes a 
single fiare event, in which one population of electrons is 
injected into the jet and causes a sudden outburst of radi- 
ation. We should also note that we treat only the optically 
thin case here. 

We will begin with the calculation of the electron number 
densities in the linear and nonlinear case in section [2l We 
will also discuss the temporal evolution of the respective 
electron densities and will outline their behavior by evalu- 
ating the upper and lower cut-offs of the spectra. We will 
show that any non- linear cooling will ultimately become lin- 
ear after some time, which requires a treatment with both 
cooling scenarios in one equation. This has been carried out 
for mono-energetic electrons by SBM, and we will extend 
their research here, as well, in section [3] 
Sections |4] |5] and |6] are devoted to the calculation of the 
total linear and nonlinear synchrotron fluence and the flu- 
ence of the combined cooling scenario, respectively. We will 
discuss our results in section jT] 

2. Linear and non-linear electron number densities 

In this section we calculate the distribution function of the 
cooled electrons for linear synchrotron and nonlinear SSC 
cooling, respectively, and discuss their temporal evolution. 
The electron number density 71(7, t) is governed by the com- 
petition between time-varying energy losses and the injec- 
tion of relativistic electrons into the jet. This competition 
is described by the partial differential equation (Kardashev 
(1962)) 



dn{'-f,t) d 



dt 



dj 



]7|n(7,t)]=Q(7,i) 



(2) 



(5(7, t) is the injection rate and we assume an instantaneous 
injection of power-law distributed electrons in the form 



g(7, t) = 907"' H [7 - 71] H [72 - 7] S (i) 



(3) 



where s is the spectral index, 71 and 72 are the lower and 
upper cut-offs of the electron spectrum, respectively, H [x] 
denotes Heaviside's step function, and S (x) is Dirac's 6- 
distribution. 



2.1. Linear cooling 

The linear cooling is described by the linear energy loss 
term 



\i\sYN = '^:^^Bl'^ = ^07^ , 



(4) 



where Ub = B'^/Stt is the magnetic energy density of a 
magnetic field B = bG. Schlickeiser & Lerche (2008) ob- 
tained the solution for the linear cooling: 



nsYN{"t,t) = qaj ''{1-Dojty 



H 



71 



l + Dojit 



X H 



72 



1 + Do-f2t 



(5) 
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The temporal evolution (with y = AqI, and Aq defined in 
the next subsection) of this solution is shown in Figs. [l]-M 
for three different cases of s. 
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Fig. 1. Linear electron distribution for s = 1.5 at four dif- 
ferent times y for the initial values: 71 = 10^, and 72 = lO'*. 
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Fig. 3. Linear electron distribution for s = 3.5 at four dif- 
ferent times y for the initial values: 71 — 10^, and 72 = 10*. 



2.2. Nonlinear cooling 

According to RS the nonlinear SSC cooling rate is defined 

by 



lilssc = ^ol / d77 n(7,i) 




(6) 



with the constant ^o — 'ici(jTP{)RE]:/rac? . Here ar is 
the Thomson cross-section of electron scattering, R is the 
radius of the spherical source, Pq = a//2-\/3?i = 3.2 • 
IQi^eV-is-i, Ek = 1.16- 10-^6 eV, ci w 0.684, and mc^ is 
the rest-mass energy of an electron. 

The integral over the electron number density 71,(7, ^) ba- 
sically introduces a time dependence of the cooling rate. 
This is reasonable, since cooler electrons have less energy 
they can transfer to the photons, which implies that the 
electrons do not cool effectively any more. 
The solution to the differential equation (pi is derived in 
appendix [X] and becomes 



"SSC (7, r) =9o7 'O'-iTY 
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72 



X H 



7 
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(7) 



The implicit time variable T is defined by the nonlinear 
differential equation 



Fig. 2. Linear electron distribution for s = 2.5 at four dif- 
ferent times y for the initial values: 71 = 10^, and 72 = 10*. 
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As one can see, the electron distributions piles up at the 
high end of the spectrum for spectral indices s < 2, as it 
was already discovered by Kardashev (1962). This does not 
happen for larger s, and in fact these cases show the exact 
opposite behavior with a drop at the high-energy end. It 
is also obvious that s only determines the steepness of the 
spectrum, and not the cooling time scales. 



where y — AqI. As we will see later the derivative of T with 
respect to y is as important as the complete time variable 
T itself in order to calculate the synchrotron spectra (cf. 
section [5]). The calculation of the former is presented in 
appendix |B1 and becomes for s > 3 



Uix) 




s-3 92 

s— 1 x^ 



{s-l)x^ 



[i-^n 



< a; < g^ 
92^ <x<l 
x> 1 



(9) 
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s-3 92 
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s — 1 x^ 



9\- 



< X < 32-1 
32-1 < a; < 1 
x> 1 



(10) 



- (s-l)2;2 

with the new time variable x = jiT, the ratio between the 
initial cut-offs 32 = 72/71, and the constants Uq 



gall 



s-3 



and Ui 



90 7i 



2.2.1. Temporal evolution of the nonlinear electron number 
density 

If one is interested in the actual temporal evolution of the 
electron distribution, one has to perform the integration of 
U{x) with respect to the time-variable, yielding for s > 3 



liUo [1 

iiUo[y- 



-5n 

'3/ 



37if/ofEf(l- 
and for 1 < s < 3 



92-liy-c3) 



1/3 



71^:^1 [92 



3-s 



1 y 



27it/i!5f(y-c5) 



37it/if5f(l 



.92 liy-ce) 



1/3 



x{y) = 

0<y<yi 
yi<y <y2 

2/ > 2/2 

x{y) = 

0<y<y3 
y3<y <y4 . 
a; > 2/4 



(11) 



(12) 



The integration and the approximations used for it can be 
found in appendix O The time limits yi and the constants 
Ci are chosen in such a way that the solution x{y) is con- 
tinuous for all y. The values are also listed in appendix [C| 
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Fig. 5. Nonlinear electron distribution for s ~ 2.5 at four 
different times y for the initial values: 71 = 10^, and 72 — 
10*. 
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Fig. 6. Nonlinear electron distribution for s = 3.5 at four 
different times y for the initial values: 71 = 10^. and 72 — 
10". 



electron spectrum for s = 1.5, s = 2.5, and s = 3.5, respec- 
tively. 

As one can see, the spectrum cools significantly faster for 
harder spectra (smaller spectral index) . It is also quite obvi- 
ous that the spectrum contracts rapidly, which means that 
for late times the spectrum is identical to a (5-function (cf. 
next subsection) . One should also notice the increasing pile- 
up for s < 2 at the upper limit for later times, which be- 
comes a drop for larger s, similarly to the linear case. 



Fig. 4. Nonlinear electron distribution for s = 1.5 at four 2.3. Temporal evolution of the cut-offs 

different times y for the initial values: 71 = lO^, and 72 = xhe qualitative behavior outlined in the last subsections of 

the electron densities is confirmed by Fig. |7] showing the 
time dependence of the upper and lower cut-offs defined by 
the Heaviside functions of the respective electron densities. 



In Figs. |4]-[6]we plotted the time evolution of the nonlinear Since the nonlinear cooling depends critically on s, we plot- 
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Fig. 7. Time dependence of the upper (72, upper curves) Fig. 8. The same as in Fig. [t] but with i?i5 = 10, and 



and lower (71 , lower curves) cut-offs of the electron spectra, 
respectively, for the standard values 71 — 10^, 72 = 10^, 
i?i5 = 1, and (75 = 1. Full curve: Linear electron cooling for 
s = 2. Dashed curve: Nonlinear electron cooling for s = 2. 
Dotted curve: Nonlinear electron cooling for s — 3.5. 



(75 = 10 . Note the different values on the horizontal axis. 



the number density of the combined linear and non-linear 
cooled (hereafter referred to as "combined" cooling) elec- 
trons in the next section. 



ted it for the two cases s = 2, and s — 3.5; for the linear case 
we used s = 2, but as we have seen before this choice is not 
important. For early times the cut-offs of the nonlinear sce- 
nario remain practically unchanged. However, the cooling 
of the high-energetic electrons sets in about three orders of 
magnitude earlier than the cooling of the low-energetic ones 
resulting in the quenching of the spectrum, which leads to 
an electron spectrum at later times that basically looks like 
a ^-function. This is expected from the denominator of the 
Heaviside functions in the electron number density, because 
of the factor 172 — 10"^ in these examples. It is also obvious 
that at these later times the nonlinear electron distribu- 
tions cool with a j/~^/^-dependence. In the intermediate 
time range the cooling of the softer spectrum (larger spec- 
tral index) is faster than for the harder case. 
Concerning the linear cooling case, we find that once 
the cooling has begun, the electrons cool with a y~^- 
dependence as expected. Similarly to the nonlinear case, 
the cooling of the high-energetic electrons begins a factor of 
about g2^ earlier than for the low-energetic ones, resulting 
again in a (5-function appearance of the electron spectrum 
for later times. 

In Fig. [7] the linear cooling operates much faster than the 
nonlinear cooling. This depends critically on the param- 
eters we have chosen for the plot: R = 10^^i?i5 cm, and 
go = 10^95 cm~^. Especially for higher values of these pa- 
rameters (cf. Fig. Is]) the linear cooling begins at much later 
times for larger R, or the nonlinear cooling begins much 
earlier for larger qo (the magnetic field strength does not 
change the overall result, since both cooling scenarios de- 
pend on B^; cf. the definitions of Dq and Aq). We note, 
however, that at one point in time the linear cooling will 
always dominate over the nonlinear cooling, since it cools 
faster for later times than the nonlinear case. 
This behavior indicates that linear and nonlinear cooling 
should not be treated separately, as we did in this section. 
This fact has been treated by SBM for the illustrative case 
of mono-energetic electrons (at all times) . We will calculate 



3. Combined cooled electron number density 

The combined cooling is described by the energy loss term 

lilcoM = |7|syjv + ItIssc 
^Aoj' \Ko+ I d-/j'ncoMh,t) I , (13) 



with Ko = Do /An. 

Inserting this into Eq. (p| yields the solution 



X H 



72 



7- 



1 + 727" 

71 



I-I-71T 



(14) 



The similarity of Eq. (14) with Eq. (l7| is not surprising, 
since the way to derive the solution is the same in both 
ways. The important difference is that the implicit time 
variable T has changed to the new definition r: 



dr 
dy 



Kq+ d7 7^ ricoM (7; t") 



(15) 



In appendix [D] we derive the complete solution of the im- 
plicit time variable x — 71 r, which becomes for s > 3 



5(2/, "0 < 1) = 

ji{Ko + Uo)y ,0<y<y2 
71^0(2/ -C4) ,y>y2 



ji{Ko + UQ)y 

[SaljiKoiy- 
liKo{y-ce) 



,1/3 



x{y,ao > 1) : 

,0<y<y2 

,y2<y<y3 
,y>y3 



(16) 



(17) 
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and for 1 < s < 3 



(s-l)any-a2J 

71^0(2/ -d5) 



3ao7iisro(2/-d6)] 
^ liKoiy-dr) 



x{y,aQ < 1) = 
,0<y<yl 

,y*i<y<y2 , 
,y>y2 

x{y,ao > 1) : 
,0<y<y*i 

1 

,y*i <y <y*3 



synchrotron flux compared to the Compton flux leading to 
the hnear synchrotron cooling of the electrons, 



as one can 



see in Eqs. (16 1 and (18). Similarly, an injection parameter 
being larger than unity leads to the non-linear SSC cooling 



(18) 



1/3 



(19) 



yl<y <yl 
y>yl 



The approximations and the values of the constants yi , y* , 
Ci , and di can be found in appendix ID] as well. They are 
chosen in such a way that x is continuous for all times 
There we also introduce (3 — (3 — s)(l 
important injection parameter 



gl^")/2 and the 



shown in Eqs. ( 17 1 and ( 19 1 , since the flux of the Compton 
peak is higher than of the synchrotron peak, which was the 
condition for the non-linear cooling. One should note, how- 
ever, that our remarks from the previous section hold: Any 
non-linear cooling will sooner or later turn into linear cool- 
ing. 

It is also noteworthy that ao is independent of the mag- 
netic field strength. This results from the same dependence 
of both, Aq and Dq, on b^. 

3.2. Temporal evolution of the combined electron number 
density 

As we did in the separated cases we also present here the 
temporal evolution of ncoAi- We will begin with the case 
ao ^ 1 and then proceed with the opposite case. 



ao 



'9071 



'"^(1-32^^) 



{s - l)Ko 



(s - l)Ko 



(20) 



for extended power law injection spectra. The power law in- 
jection parameter ag depends on the particle number den- 
sity go I the initial lower cut-off 71, the electron spectral 
index s, the constant Kq, and only weakly on the upper 
cut-off, since gl~'' <C 1. 

3.1. Power law injection parameter 

For standard blazar parameters qo = lO^gs cm~'^ and R = 
10i5i?i5 cm we use from SBM (Ko/qo)'^/^ = lOGiRi^qs)-^/^ 
and introduce the total number of instantaneously injected 
electrons 



N 



10^°N5, 



An 



-R\o 



^77 
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•^- 



~-~s 
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(21) 



yielding q^ = {s - l)N5o/iA.2Rl^jl- 
law injection parameter becomes 



Fig. 9. Temporal evolution of the combined electron num- 
ber density for s = 1.5 and ao = 0.1 at four different times 
"*), so that the power y with the initial values 71 = 10 and 72 = W. 
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-3 7lA^50 

-R15 



(22) 



in terms of the characteristic Lorentz factor introduced by 
SBM 



7s 



217i?i5 



N, 



1/2 



(23) 



The injection parameter is determined by the total number 
of injected electrons and the size of the source as indicated, 
but independent of the magnetic field strength. If the lower 
cutoff of the injected power law 71 is higher (smaller) than 
7b, the injection parameter ag will be larger (smaller) than 
unity. For a compact source with a large number of injected 
electrons, the injection parameter is much larger than unity. 
Basically the injection parameter is the fundamental order- 
ing parameter that defines the ratio between the magni- 
tudes of the synchrotron and the Compton peak. An injec- 
tion parameter smaller than unity corresponds to a higher 



The results are indeed comparable to the results of section 
|2.1[ The width of the electron distribution shrinks linearly 
and similarly for all shown cases. The major difference is the 
pile-up for s < 2 at the high-energy end of the distribution 
and that the spectra steepen for higher values of s. 



The important point of the plots [12] -[T4| is that the cooling 
is much quicker than in the previous case for ao <C 1. One 
should also notice the pile-up again at the high-energy end 
of the electron distribution for the hard spectrum, while for 
higher spectral indices the spectrum drops off. 

3.3. Temporal evolution of the cut-offs 

In this section we want to verify the results of the last sec- 
tion, especially of the temporal behavior of the electron 
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Fig. 10. Temporal evolution of the combined electron num- 
ber density for s = 2.5 and ao = 0.1 at four different times 
y with the initial values 71 = 10 and 72 = 10^. 
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Fig. 12. Temporal evolution of the combined electron num- 
ber density for s = 1.5 and ag = 100 at four different times 
y with the initial values 71 = 10 and 72 = 10*. 
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Fig. 11. Temporal evolution of the combined electron num- 
ber density for 3 — i.b and ao = 0.1 at four different times 
y with the initial values 71 = 10 and 72 — 10*. 
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Fig. 13. Temporal evolution of the combined electron num- 
ber density for s = 2.5 and a^ — 100 at four different times 
y with the initial values 71 — 10 and 72 = 10*. 



distribution. Therefore, we discuss the evolution of the cut- 
offs for three different cases of s shown in the Figs. [Tsjand 

m 

Since for ao *^ 1 the behavior of the cut-offs is independent 
of s, we only plotted one case in Fig. [15] It is obvious that 
this is indeed very similar to the linear solution of section 
(with the plot of the cut-offs in section 2.3). 



2.1 



For ao 3> 1 we notice that the cooling begins much ear- 
lier than in the other case, with the hard case beginning 
even two orders of magnitude earlier in the example shown 
here than the soft case. We also notice that the soft case 
cools linearly in the beginning until it merges again with 
the already non-linearly cooled hard case (here, at roughly 
y K, 10""'^*). Interestingly, form that moment on both cases 
cool the same. When both begin to cool linearly again at 
late times (at about y > 10^^) they match the cooling of 
the case with ao <C 1. 



Obviously, the electron distributions of all cases undergo se- 
vere quenching resulting in the appearance of a (5-function 
for later times. 

Overall we can say that the behavior of the combined cool- 
ing is similar to the respective cases of the separated cooling 
scenarios. In order to stress the most important point once 
more, the non-linear case turns linear at late times as we 



expected (cf. section 2.3) 



4. Linear synchrotron fluence 

In the next three sections we will derive the synchrotron 
fluence spectra, for which we need at first the respective in- 
tensity spectra. We use the monochromatic approximation 
of the synchrotron power introduced by Felten & Morrison 
(1966). In general, the synchrotron intensity spectrum of a 
homogeneous, spherical, optically thin source is then given 
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Fig. 14. Temporal evolution of the combined electron num- 
ber density for s = 3.5 and ao = 100 at four different times 
y with the initial values 71 = 10 and 72 — 10^. 




log{y) 

Fig. 16. Temporal evolution of the cut-offs of the combined 
electron number density for s — 2.0 (full curve), s = 3.5 
(dashed curve) and ao — 100 for the initial values 71 = 10 
and 72 = 10"*. 



time variable x — jit — -j-y, we obtain 



IsYN{f,x) 



RgocaxB'^ 



X H 



./- 



1 



(l + ^ox)^ 



H 



1 



(32^' + Doxy 



-f 



(26) 



The fluence is the time-integrated intensity spectrum. 
Setting an arbitrary upper limit would result in the frac- 
tional fluence, which shows the fluence after some time. 
Here we treat, however, only the total fluence with an up- 
per limit that is infinite. Thus, 



-24 -22 -20 -1 



log{y) 



F. 



SYN 



if) 



RgocaxB 



487r2i.,7f 



^/^/dx(l-i?o/^X 
i^slt J V 



s-2 



Fig. 15. Temporal evolution of the cut-offs of the combined 
electron number density for s = 2.0 and ao — 0.1 for the 
initial values 71 = 10 and 72 = 10'*. 



X H 



/- 



1 



(l + Dox)' 



H 



1 



ifh+Doxy 



(27) 



Introducing the variable fi = DqJ^^^X "^^ yield 



by 



/(z.,i) = 



R caxB'^ 
47r Gtt 



F< 



SYN 



(/) = 



d7 n{j, t)^"^ 5 [v — ^sl"^) 



(24) 



RqocaxB'^ 
487r2z/37f/«/2_D, 



00 

- / df7(l-r2)'^ 
'0 J 



X H 



f^- (1-/1/2)1 H[(l-.9XfX-f^ 



(28) 



where v is the frequency of the emitted photon. The char- 
acteristic frequency is defined by Vg — ..ffj^^a — 4.2 -10^6 Hz. 



47r7nc^ 



We can use the substitution rule for 5- functions, providing 
271/3 (5 (i/ — 1/57^) = 5(7— yjvjvs ) . Inserting this into 
Eq. (p4| results in: 



The result of this integral gives two different frequency 
regimes, and can be written as 



I{v,t) 



R cotB'^ 
in 127rz/, 



d7n(7,i)7 S 7 



(25) 



If we now apply this to the linear electron distribution Eq. 
5|, define the normalized frequency / — — ^, and a new 



f-"''{i-sl-) 



with Ni = ■„ f''"?;^^ ,, . 



FsynU) = Ni 



(29) 



This is the well known result that has been obtained by 
Schlickeiser & Lerche (2008), too. The important points of 
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this result are that the fluence is independent of the spec- the ongoing discussion into two parts for each case of the 
tral index s of the electron injection rate for frequencies spectral index, 
smaller than unity following a power-law with a spectral 



index asYN = 1/2- At / = 1 the spectrum exhibits a break 
to a steeper spectrum with Aa syn — ^-j^ depending on s. 
The spectrum is cut off at / = ^f. 



5.1. Large spectral index 

We begin with the spectral index being larger than 3. Then, 
the first integral we have to calculate is for the case < 
UJ < g2^f^'\ where U{w/f^/^) = Uo{l - gt'), yielding 




Fsscif) = A^o 



-91- 



9-.'l"' 



da;(l-w)" 



X H 



(1-/V2)1 llhl-g-^f'^) 



(32) 



where we set iVo = /f'^^^ . 

The integration can be easily performed giving, according 
to the Heaviside functions, solutions for two regimes of the 
normalized frequency: 



(i) 



^sscU) = No 

s-l 



f-s/2 



1 



l-9t'-^-^ 



1/2 



Fig. 17. NFgYN as a function of / for three cases of s {Full: 
s = 2, Dashed: s = 3, Dot-dashed: 3 — 4) and for 02 = lO''. 



1 < / < 35I 



5i</<.9i 



(33) 



The red hues indicate the analytical solution Eq. (29 1 with The intermediate w-range (g^^ f ^/^ < ^^ < fi/S) jg the next 
an offset of 10-0-7. B yy2 J - - J J 



case we have to consider. Here 

2(a.//V2)-3 



1- 



l-s 
92 



1 



s-3 

s - 1 (a;//i/2)2 



(34) 



In Fig. [17] we show the numerical solution of the flu- U{uj) = Uq 

ence spectrum for which we normalized the fluence by 

NFsYN = Fsyn/Ni. We plotted the result for 5-2 = 10^ Approximating this with the leading term, which is valid 

and three values of s. We also plotted the analytical re- foj. j^^^gt p^rts of the time period and for spectral indices 

suits, but with a small offset in order to highlight them. ^ot too close to 3 (cf. Fig. ml in appendix O, the integral 

As one can see, the numerical and the analytical result becomes 

match each other rather well, which is not so surprising, 
since we introduced no strong approximations during the 
analytical work. 



/ 



1/2 



5. Nonlinear synchrotron fluence 

In order to obtain the intensity spectrum and the fluence of 
the nonlinear case we have to insert the nonlinear electron 
distribution Eq. Q into Eq. ( 25 I , yielding 

W(/,.) = ^r?^7r/-(l-/^/^.^^"^ 



Fsscif) ^ Nof-'/^ J dLuil-Luy-^ 

9^'r^' 
xhL-(i-/V2)1 uUi-g-^f/^)- 



(35) 



This is again easily solved for two regimes of / obtained 
from the Heaviside functions and indicated below. We find 



AStt'^i^, 



Fsscif) = Nof 



-s/2 



1 



X H 



f- 



1 



(l + a;)2 



H 



1 



(32 



-/ 



(30) 



Integrating this with respect to time yields the total fluence: 

Fsscif)^ :7'r\" j -'^' I duju-^ 




s-l 
, 1/4 < / < 1 

,l</<3.9i 



(36) 



x{l-uj) 



s-2 



H 



Jl/2 
'^-il-f'") 



The third time regime lo > /^/^ requires U = C/ofufll 
g\~'')fuj~'^ , for which the fluence becomes 



Fsscif) - iVo^(i - gl-T'f-"^ 

s — 6 



dww^ 



X H 



il-92'f''') 



(31) 



for which we used the substitution lo = f^l'^x. 

Since we must take care of the spectral index s, we split 



x(l-w) 



X H 



s-2 



H 



^1/2 
(1-/^/^) 



(1-52-V^/^) 



U) 



(37) 
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The primitive of the integral can be obtained by two inte- 
gration by parts, yielding for arbitrary limits a and b, which 
need to be specified by the Heaviside functions later, 



-s/2 







-\b 



S-1 



(1-^) 



s-1 



2io 



■{l-uY 



-(l-c) 



s+l 



(38) 



For / < 1/4 the limits become a — 1 — /^/^ and b — 
(1 — (72^^/^/^). Assuming 52 ^ 1 we can neglect the negative 
terms, which are the terms where we inserted the upper 
limit. Since / < 1/4 we see that the second and third of 
the remaining terms are also much reduced compared to 
the first one. Thus, we approximate the integral with 



INT w ^ '-f- 

s — 1 



(39) 



From the Heaviside functions of Eq. (37) we obtain another 
frequency regime located between l/T^J / < 1- The limits 
in this case become a = /^/^ and b = l — g^^f^^^. Using the 
same approximation as above, for which we note however 
that these approximations are not as well fitting as in the 



Collecting terms we find for the total fiuence with a large 
spectral index s > 3 



FsscU) - N^ [Fiif < 1/4) + F2(l/4 < / < 1) 
+^3(1 < / < ig2) + F.i^gl <f< gl) 



with the terms 
Fi 

j-s/2 , 



[l-f'^f .^3/2 

{s-i){l~gl~^y 



F 



s-1 



/^_(l_/l/2)-l 
y-s/2 



(l-/l/2)s-l 

{s-3)il-gl-n' 



(42) 



(43) 



(44) 



/- 


5/2 


(,s-l)(l-.g3-) 


f-s/2 
<-l 


Yi-(^ 

V V32 



1- 1 



1/2N 



1/2N 



s-1 



Fi 



1- 

j-s/2 



{s-l){l-g-l 



3-SN 



(45) 
(46) 



In the third part we approximated again for 52 3> 1. If 
we do the same in the fourth part, the third and fourth 
part would be the same and could be merged. This does 
not change much to the overall result, which should be dis- 
cussed briefly. 

The fluence behaves like a power-law with a constant spec- 
tral index assc = 3/2 for / < 1, independently of s, which 
is steeper than in the case of linear cooling. At / sa 1 the 
spectrum exhibits a spectral break to a steeper power-law 
with a change in the spectral index of lS.assc = {s — 3)/2. 
At a frequency / = g| the spectrum is cut off at the high 
energy end. 



iNT« ^ (1 f/^y ' . 


(40) 


A 


Thus, we yield in the late time case for the fluence 




LL 

2 


Fsscif) = No'-^a - gi-T'r"^ 

s - 






/ /I rl/2N2 =_i 

1^(1-/^/^)^" ,l/4</<l ■ 


(41) 






Fig. 18. NFssc as a function of / for three cases of s > 3 
{Full: s = 3.1, Dashed: s = 4, Dot-dashed: s — 5) and for 
02 = 10^. The red lines indicate the analytical solution Eqs. 
(43l - (lief with an offset of 10""-^. 



In Fig. [18] we show the numerical solution of the flu- 
ence spectrum for which we normalized the fluence by 
NFssc = Fssc/Nq. We plotted the result for 52 = 10^ 
and three values of s. We also plotted the analytical re- 
sults, but with a small offset in order to highlight them. 
Obviously for most parts of / the analytical result and the 
numerical result match each other. However, for 1/4 < / < 
1 there is a discrepancy between the results, which is also 
quite obvious, since the analytical result is not continuous 
at / = 1/4. We already indicated during the analytical cal- 
culations that the approximation for the second part (Eq. 
( [44| ) could be invalid. One can also see that the approxi- 
mation in the third part (Eq. (45l) becomes less valid for 
s — > 3. Thus, the numerical result conflrms the expecta- 
tions we mentioned during the derivation of the analytical 
result. 
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5.2. Small spectral index 

Now, we want to deal with the fluence for electron spectral 

indices 1 < s < 3. As a matter of fact, the steps are quite 

similar to the case before, but nonetheless we will repeat 

them here. 

We begin again with the case < cj < g2^f^^^, where 

C/(c^//V2) 



Uiigt^^l): 



f-s/2 

FsscU) ^ N^^^j— 



92 f 



lfl/2 



52 



1 



duj{l--ujY 



X H 



-{l-f'H Yi\{l-g^'f/')- 



(47) 



with the definition iVi = -R|oco-t-B Since the integral is 
easily obtained once more, for small u> the fluence becomes 

j-s/2 1 



FsscU) = Ni 



91- 



ls-1 



1 



(^) 



1/2 



fe) 



1 < / < hi 



, hi <f<9l 



(48) 



where the Heaviside functions of Eq. ( 47 ) defined the fre- 
quency regimes. 

The intermediate w-range (g^^/^^^ < i^ < /^^^) is the next 
case we turn our attention to. Here 



U{uj) = Ui 



2{L^/f 



l/2\s-3 



- 1 



9'2-' 



3-s 

^^(c.//l/2)2 



(49) 



Taking into account only the leading term of U{uj) (the 
validity of this approximation is shown in Fig. |C.2| in ap- 
pendix O, as we did in the same case for the steep electron 
spectrum, we obtain 



F. 



SSC 



if) = N,r''^ 



V' 



dujuj^-^il-uY 



g^'P'^ 



X H 



(1-/V^) H (1-52-7^/^) 



^-lfl/2^ 



, (50) 



This integral can be expressed in terms of the hypergeomet- 
ric function, yielding no simple solution. Thus, for interme- 
diate frequencies we cannot present an analytical solution. 
Finally, we have to turn our attention to the late time 
regime, which means uj > f^^^. Here U — C^ifiif ^(1 ^ 
gl~''), for which the fiuence becomes 



Fsscif) = Ni 



f~ 



1 — S Q 

.92 -^ 



dujuj (1 — Lj) 



s-2 



fl/2 



X H 



■..-(1-/1/2)1 h[(1-52-V^/')-^ 



(51) 



This is almost the same integral that we solved in the case 
of s > 3, thus we can use the same primitive (Eq. (38 1). 



In fact, the limits and the frequency regimes are the same 
as in the case of steeply injected electrons. Therefore, using 



the same approximations as we did in that case, we obtain 
the solution for / < 1/4 



INT « -^ -f — 

s — 1 

and for 1/4 < / < 1 



INT: 



/ 



1 



1-/ 



1/2 



(52) 



(53) 



Collecting terms, we find the fiuence for large oj: 



Fsscif) = Ni 
-/' 



^(1-/ 



1/2 



s-1 



:ii- 9l-T'f-^ 

, / < 1/4 
, 1/4 < / < 1 



(54) 



Now, we have performed all necessary steps to calculate 
the synchrotron fiuence of SSC cooled electrons that are 
injected into the jet with a spectral index 1 < s < 3. 
Therefore, 



Fsscif) =Ni 



with 



F5if<l/^)+Fei^gl<f<gl) 



(1 - /V2)2 

i'F. - T^irJ 



F« 



(3-,s)(l-.9r^) 

/-/^ 
(,5-l)(.g3--l) 



l-'4 



(55) 

(56) 
(57) 



We do not present the results for the intermediate frequency 
regime 1/4 </ < 1/4(7|, because this is the part where the 
intermediate cj-regime is also important, if not more impor- 
tant than the contributions of the early and late w-regimes. 
And for that part we did not obtain an analytical solution. 
However, the results in the other frequency parts are note- 
worthy. The fiuence for fiat injected electron spectra shows 
similar characteristics as the fiuence for steep injected elec- 
trons. We find, again, that the spectral index does not de- 
pend on s for frequencies / < 1/4 and has the same value 
otssc = 3/2. There is also a cut-off for high frequencies at 

f = 9l 

In Fig. [19] we show the numerical solution of the flu- 
ence spectrum for which we normalized the fluence by 
NFssc = Fssc/Ni. We plotted the result for 52 = 10^ 
and three values of s. We also plotted the analytical re- 
sults, but with a small offset in order to highlight them. 
As one can see, the numerical and the analytical solu- 
tions agree very well for the presented frequency intervals. 
Interestingly, there is no spectral break in the numerical re- 
sult at around f — I, different from the behavior of steeply 
injected electrons. The spectrum is cut off at / = 5I as 
expected. 

It is interesting that the spectrum does not exhibit a break 
at / = 1, implying that the spectrum for / > 1 does not 
depend strongly on s, either. It is quite obvious that the 
offset between the numerical and the analytical solution 
does not change from the low frequency regime to the high 
frequency regime near the cut-off. The curves of the ana- 
lytical plot at the cut-off are also equal for all cases, similar 
to the numerical result. Especially this last point is differ- 
ent from the steep injection case, where the magnitudes of 
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6.1. Large spectral index, small injection parameter 

Here we discuss the case of s > 3 and ao ^ 1- 
The first regime we have to consider is for < ^t < ^2, 
with *2 := *(j;2) = /^/^. With d* = P^^jdKo + Uq) dy 
we find 



FcOmU) — Fa 



f-s/2 



d*(l-*)" 



X H 




(1-/'/^)-*] H [vl- - (1 - /I/2; 



Fa 



ji{Ko + Uo){s- 



lyha-^M^ 



(60) 



12 3 

log(f) 

Fig. 19. NFssc as a function of / for three cases of s > 3 
{Full: s = 3.1, Dashed: s = 4, Dot-dashed: s — 5) and for 
Q2 = 10^. The red lines indicate the analytical solution Eqs. 



where a and b are defined by the limits of the integral and 
the Heaviside functions. For / < 1 we find a — 1 — /^/^ and 
b — /^", resulting in 



FcOMif<l)^Fa 



f 



-1/2 



7i(i^o + t/o)(s-l) 



(56) and (57 1 with an offset of 10 



I -if 



-1/2 



\s-l 



(61) 



. . This becomes negative for / < 1/4 which is not physically 

the curves differ significantly from case to case. These are g^jble, and, therefore, this solution does not contribute 

hints tor the validity of the numerical solution exhibiting r- i r 

lor lower ireQiieiicies 

f„l^iL^^'„'^„^i''^ '" ""'' important difference from the steep p^^ ^ > ^ ^^ obtain 'the limits a - and 6 - 1 - p^'g^\ 



injection case. 

6. Synchrotron fluence of combined cooled 
electrons 

We will; now, calculate the synchrotron fluence spectrum 
of electrons undergoing the combined cooling. 



resulting in 

FcOAlif > 1) = Fa 



f-s/2 



7i(i^o + f/o)(s-l) 
/ 



As before, we begin with the intensity by inserting Eq. (|14|) 
intoEq. (|25|. With th 
frequency / we obtain 

qoRcarB^ 



1 - 



(62) 



into Eq. ( p5| . With the same deflnition for the normahzed Obviously, this is cut off at / = ^f. 

For * > v]/2 we flnd d* = f^^^-fiKo dy , leading to 



Icom{Lx) = 
X H 



1 



Wrf'^i}-!"'^) 



-./ 



H 



FcOAlif) = Fa^-^ I d*(l-*)'^-2 



./- 



(l + i) 



^^2 



liKo 



(58) 



X H 



The integration with respect to time yields the synchrotron 
fluence: 



f-s/2 



^Fa 



liKois - 1) 



H *- (1-/1/2) 



(63) 



FcoAiif) = / dtIcoM{f,S:) = j- dy IcoM{f,x{y)) 



00 

^„/^ / dy (1 - f^'^i{y)y'^ Fcouif < 1/4) = Fa 



We flnd for / < 1/4 that a = 1 - /^/^ and 6 = 1- f^^^g^^ 
yielding 



= Fa 



I 



-1/2 



ly [1-52-1 



X H 



ig2^ + i{y)? 



-.f 



H 



/- 



1 



m 



with the definition _F^ 



go Rc(TtB 



liKois-l) 
For / > 1/4 we obtain for a = /^/^ and 6=1- f^^'^g2^ 

FcoAiif > 1/4) ^ Fa ^ 
X "(r V2 _ 1 



(64) 



71^0(5-1) 

s-l 



.9^^ 



It is advantageous in this case to use y as the integration 

variable instead of a;, as we did in the other cases, because 

we already made a lot of approximations while calculating which becomes negative for / > 1. 

x{y) for the discussions in section [3] and it is unnecessary Collecting terms, we obtain the total fluence 

to repeat these steps again. \ it 

Using i{y) form Eqs. ([16]) - ([l9| we will in each of the fol- Fcom{s > 3, ao < 1) = F„ [F^ (/ < 1/4) 

lowing integrations deflne the new variable ^ := f^/^x{y). +F^(l/4 < / < 1) + F^{1 < f < (72)] ■ 



(65) 



(66) 
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with 
















F^ - 


j,Ko{s 1) ^' ^2 J 
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/-1/2 ^ 1 


"(/-^/^ - 
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-9\-'' 
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■l_(/-i/2_ 


_^y-l 
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F! = 


7i(A'o + C/o 
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-1) 


1- 




s-1 n 
2 





This is, again, cut off at / = gf. 

For ^F]; < <!> < **, where vp^ = v[/(y*) = /i/2^ the differen- 
(67) tial becomes d* = /^/^ r^l'n ^2 dy , leading to 



J^COA/Cn^J-o ^'t^' J^"° / d*(l-Vl/)-2 



(68) 



(69) 



This is the expected result. Since the cooling is dominated 



l^iUiP 



*r 



X H 



(1 - /i/^') - *1 H [* - (1 - /1/2) 



F, 



271 [/i/3 



(1-*)^ 



(74) 



by the linear part, the solution is comparable to the linear y^ n -c • j- ^ 1 j-i. r -j- u 1 -fi/2 

„■' r • i-n • 11 1 . , i'or small irequencies f < 1 the limits become a = 1 — / ' 

nuence 01 section |4| representing a broken power-law with Ah— fi/2 It 

a constant spectral index {j3i — 1/2) for frequencies below ' ' 



unity, and a steepening of A/3i = (s — l)/2 for / > 1. The 
spectrum is cut off at / = 5I . 

6.2. Small spectral index, small injection parameter 

We continue with the case of 1 < s < 3 and ap ^ 1- 
The first part is < ^ < ^^, with ^l := ^{yl) = /^/^^^^ 
In this case, we find d* = f^/'^^i{KQ + Uigl~'')Ay leading 
to 



FcomU <l)^F^ 



/ 



-1/2^,2 



271C/1/3 



(/- 



/2 



(75) 



This turns negative for / < 1/4 

In the case / > 1 we obtain a — /^^^g^^ and h — \ 

/1/25-1. Therefore, 



FcomU >l)^F^ 
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^1/2^2 



a, 







FcomU) ~ Fa 



f-./2 



*r 



271C/1/3 



7i(ifo + t/i5r') 



d*(l-*)" 



1- 



1/2N 



(76) 



X H 

Fr 



(l-/'/V)-*J H [vl/ _ (1 _ /V2) 

^^ [-(1 - *)^-l]' 

7i(i^o + c/i.gr^)(s-i)^ ^ ^ ^-^ 



which does also not contribute for / > 1/4(;|, since 
FcoM < in that part. 
(70) For * > *5 we find d* = P^'^JiKq dy , leading to 



With the help of the integration limits and the Heaviside 
functions we determine a and b for three different regimes 
of /. For / < 1 a = 1 - /^/^ and b = f^/^g^^ for which the 
fluence becomes 



FcomU) = F, 



f 



-s/2 



liKa 



d* (1-^^-2 



*? 



FcomU < I) ^ F^ 



f 



-1/2 



X H 



(1 - f^^92') - *1 H [* _ (1 _ /1/2) 



7i(ifo + C/i52"^)(s-l) 



l~f- 



(71) 



F, 



r 



-s/2 



"71^:0(5- 



lyHi-^M: 



(77) 



Since / < 1 and s > 1, this is always negative, which The results are obvious^^ similar to the st^ep inject wii case, 
means, in turn, that this part does not contribute to the s° that for / < 1/4 a = 1 - / and 6=1-/ ^2 t^e 



overall result. 

For 1 < / < 1/4(7| we obtain the limits a ~ Q and b 

f^''^g2^, resulting in 

FcoMil <f< 1/452') = Fa 



result becomes 



FcomU < 1/4) = ^o .^ ,; ,„ ^. [l - 52^1 



7i(i^o + C/iff2'"')(s-l) 
-1' 



'71^0(5-1) 



(78) 



2„-l 



1 



1/2N 



For / > 1/4 we obtain for a = /^/^ and 6=1- ,f^^'g2 
(72) FcoM U > 1/4) = Fa —f- 



The third interval is for / > l/4g| with a — and b 
1-/1/25-1. Thus, 

i^coM(/>l/45|)=^a- ^ ' 



-1/2 



(r^^^- 1)^-^-52^ 



l-s 



(79) 



which becomes negative for / > 1. 

Collecting terms, we obtain the total fluence for this case 



7i(Ao + C/,.<;r')(»-l) 



52 



(73) 



Fcom{1 < s < 3,ao « 1) = Fa [F^U < 1/4) 
+FH1/4 < / < 1) + F6(1< / < 1/4.92) 

+i^J(l/452 < / < g2)] , 



(80) 
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(81) 



(82) 



F^ = 






(84) 



In F^ we inserted the definitions of ag and /3. which made 
it possible to find F^ = F^. 

The main results are similar to the steep injection. For / < 
1 the spectrum shows a power-law with a constant spectral 
index /3i = 1/2, while for / > 1 the spectrum steepens with 
A/3i = (s - l)/2. We also find the cut-off at / = 5|. 

6.3. Numerical solution for a small injection parameter 

Before we proceed with ao S> 1, we want to compare the 
analytical result with a simple numerical integration, as we 
did for the linear and non- linear cooling, as well. For the 
plots we dropped the constant F^, which means that we 
normalized the fiuence by NFcom = Fcom/ Fa- 
in the following plots we fix ap, and, thus, q^ is not a free 
parameter any more, but becomes 



go = 



alKo{s-l) 



7^(1 



9l- 



') 



(85) 



This confirms what we discussed in section |2.3| that go 
strongly infiuences which type of cooling occurs. 

The results shown in Figs. |20] and [21] do not match each 
other as well as they did in the linear and non-linear cool- 
ing cases. Although we inserted no offset, the analytical 
solution has a smaller magnitude than the numerical re- 
sult. 

For frequencies below unity one can see, however, that the 
behavior of the analytical and numerical solutions are the 
same. Both show a power-law with a spectral index of 
/3i = 1/2. For frequencies above unity and for s > 2 the re- 
sults also agree rather well, and show the cut-off at / = g|. 
For s < 2 the results do not match each other so well 
any more. The numerical result cannot be represented as 
a single power-law, but steepens gradually until it reaches 
a small plateau (which one could call a small pile-up) at 
about / = l/4g|, and cuts off at / = gf. This behavior 




1 2 

log(f) 

Fig. 20. NFcoM as a function of / for three cases of s 

(83) [Full: s — 1.5, Dashed: s = 2.5, Dot-dashed: s = 3.5) and 

for ao = 0.1, and 32 = 10'^. The red lines indicate the 



analytical solution Eqs. (66 1 and (80) 



^0.01 




1 2 
log(f) 

Fig. 21. NFcoM as a function of / for three cases of s 
{Full: s = 1.5, Dashed: s = 2.5, Dot-dashed: s = 3.5) and 
for ao — 0.01, and 32 = 10'^. The red lines indicate the 



analytical solution Eqs. (66 1 and (80). 



is more obvious for ao closer to unity, which could be a 
hint that the approximations are not valid enough for this 
case. Interestingly, for s = 1.5 and ao = 0.1 the analytical 
solution is not continuous at / = 1/4^2) which is a little 
bit strange, since it should be continuous as one can easily 
verify. 



6.4. Large spectral index, large injection parameter 

We will continue our analytical calculation of the syn- 
chrotron fluence with the case of ao 3> 1, and begin with 
s >3. 

For < ^ < ^2 the calculation is the same as for a ^ 1, 
which means that we can use the results obtained during 
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the calculation of that case. Thus, 



Fcom{1/^ <f<l)=F, 



f 



-1/2 



The last interval is for / > 1/4, and with a = /^/^ and 
6=1- /i/2g-i ^e find 



j,{Ko + Uo){s-l) 

1 - ( r ^/' - 1 



(/- 



(86) 



and 

Fcom{1< f < gl) -- F^ 



f U-f 



1/2 



FcomU > l/^) = F„ 2 

1/2 



/- 



s-1 



1 



^-./2 



^. 



y-s/2 



a57i^o(s- 1) 

2 






""oTi-f'^ols- 1) 



1-/ 



1/2 



(91) 



1 



This solution turns negative for / > 1. 
(87) One integration remains, which is for ^ > ^3. The differ- 
ential is substituted by d$ = f^^^^fiKody . Thus, 



The next integration will be done for ^2 ^ ^ < ^3, 
where ^3 = f^^^aQ. Hence, the differential becomes d^P = 
"^^^ p/^a^^iKQ Ay , and the fluence reads 



FcomU) — Fa 



I 



-s/2 



liKa 



d1'(l-* 



,s-2 



FcOM{f)^Fa V r. I d**2(l_^).-2 



X H 



J: 

(1 - f^^92^) - *1 H [* - (1 - /1/2) 



X H 



(1 - f^'92') - *] H [v& - (1 - /1/2) 

f-s/2 



Fa 



liKois - 1) 



Yy[-(i-*)]: 



(92) 



F, 



r 



a 9 



aoliKo 



*2 



S-1 



(1-* 



2* 



,s-l 



lb 



For / < ckq ^ the limits become a — I — f^^^ and b — 
1 — ,f^^^g2^, and, therefore. 



2 -(1-*)*-^ 



(1-*)" 



" 2 



ag 71/4:0 



*^ 



s- 1 



(l-*)'^ 



FcoAiif <a„^)=Fa 

In the case / > a^^ 
resulting in 

FcoMif > ao"^) ^Fa 



y-1/2 



71^0(5-1) 



n [1-52^-1 



(93) 



1-/^/% 



= /i/^ao, and 6 = 1 - P^^g^^ 

j-s/2 

(94) 



7i^o(s-l) 



1^92 



In the last step we used the leading term approximation 
again, which we have already used during the calculations 
of the non-linear fluence. We infer from the limits and the x 
Heaviside functions that the fluence of this part could con- 
tain three frequency regimes. which is always negative, and does not contribute to the 
The first one is for / < q:q" , and with a=l — /^'^ and fluence. 



b = f^^^a^ we obtain 



F, 



COM 



if < ao ') - Fc 



f^S/2 



a^JiKois - 1) 



/I J-l/2\2 r„2/r-l/2 „ \s-l 



(89) 



Inspecting this solution a little further one finds that this 
does not contribute to the overall fluence for / < (ap + l)"^- 
This means that this solution has only a very narrow regime 
of applicability. Therefore, we neglect it entirely. 
The second frequency interval is a^'^ < f < 1/4 with the F^ 
limits a = 1 — .f^^^, and & = 1 — ,f^^^g2^- Hence, 



Collecting terms the fluence for s > 3 and ao ^ 1 becomes 

Fcom{s > 3, ao » 1) = Fa [F^^f < a^^) 
+F:\ao' <f< 1/4) + F;3(1/4 < / < 1) 

+F;4(1 < / < g2)] , (95) 

with 

y-1/2 



p*i 






[^-g2 



FcoMiao^ <.f <l/4.)=F, 
f-3/2 



f 



-3/2 



F, 



*3 



al-yiKois- 1) 



fV^y 1-. 

32 



P*4 



al-fiKo{s - 1) 



7i(s-l) 
/ 






{Ko + Uo)Ko 



-s/2 



7i(ifo + C/o)(s-l) 



52 



(96) 
(97) 

(98) 
(99) 



Fa 



ao7i^o(s - 1) 



Ty[i-3n • 



What we have found is a spectrum that is represented by 
, , mainly three power-laws, instead of two as in the previous 
*- ' (especially the linear and non-linear) cases. For very low 
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frequency the spectral index is /3i = 1/2, identically to the This is identical to Eq. (88 1 and we used the same approx- 

linear cooling as expected. For intermediate frequency up to imations, again, as in that calculation. Therefore, the final 

/ « 1 the spectral index becomes /32 = 3/2, which is also results will be the same, as well, and we c an just copy them 

expected since this is identical to the non-linearly cooled from Eqs. (90 1 and (91 1, neglecting Eq. (89 1 as before: 

case. In the region around unity the spectrum changes its 

spectralindex according to Aft = (s — 3)/2. The spectrum ^ / -2 ^ r ^ t / a\ n / 

is cut off at / = gi. Fcomia, < / < 1/4) = F„ 



aljiKois- 1) 



6.5. Small spectral index, large injection parameter 

Finally, we deal with the case 1 < s < 3 and ctQ ^ 1. 
For ^ < vJ/J we can use the result obtained in section 6.2 
yielding 

FcoAiil <f< 1/452') = Fc. ^ ' 



('-/■'f-( 



92 



1/2N 



.92 



F„ 



f 



-3/2 



a^-/iKa{s- 1) 



[1 - .9M 



(105) 



and 



7l(i^o + f/l<?2'~')(s-l) 



1- 1- 



and 

^COA/(l/4ff2'</<52)=^a 



./ 



1/2N 



-s/2 



(100) 



FcoAiif > 1/4) = F„ 



a'oliKois- 1) 



/I-/ 



1/2 



f 



Fr, 



f 



-s/2 



1/2N 



f^92 



1 



"o7i^o(s-l) 



1-/ 



1/2 



(106) 



For vj/* < xp < \I/* with ^l 
tive with 



d* = *3 - sf — 
resulting in 



.(101) This solution turns negative for / > 1. 

The last integration is for $ > ^4 with the differential 

/1/2, we substitute the deriva- ^^ ^ /'^'ti^o dy . Thus, 



/ 27iL/i(4-g) 
V 5-1 



&y ,(102) 



FcoMif) = fJ-^ I d* (1 - *)-^-2 



liKo 



*,•; 



X H 



(1-/^/^52-^-* 



FcOmU) — Fa 



/^(4-s)(s-l)J^ 



d** 



s-3 



F„ 



/ 



H 

-,s/2 



*-(l-/ 



1/2N 



(s-3)(27iC/i(4-s)) 
x(l - *)-'-2 jj [(^ _ /^/ V) - * 
X H [* - (1 - /^/2 



7ii^o(s-l) 



-(l-*)]a 



(107) 






We can use the solution (93 1 from the case s > 3, yielding 

f-1/2 



(103) 



FcoMif<ao^)^F, 



liKois-l) 



IT ['-»"'] 



(108) 



The integral can be expressed in terms of the hypergeo- Since the result of Eq. (94 1 is always negative, we neglect 

metric function, but, unfortunately, one cannot obtain an it here, too. 

analytical form. Hence, we have to terminate further dis- Collecting terms we find 

cussions of this case. 

The next time step is Vl/^ < * < ^4, where ^4 w /^^^ao 

Thus, the differential becomes d^ = "^^^f^^^^aQ^iKody 

and the fluence reads 



*: 



FcOAlif) = Fa 



r 



OioliKo 



d**2(-2^_^) 



s-2 



*; 



X H 



(1 - f'^'92') - * H [* _ (1 _ /1/2) 

^2 



Fcom{1 < s < 3, ao » 1) = Fa [F^Hf < a^^) 
+F*^{a^^ <f< 1/4) + i^r (I/45I < / < gi)] , 

with 

f^3/2 



Ff 

Ff 



-f1 

~ " ,,2, 



2+s 



(^q7iKo{s 



-yy [1-52^-1 



2* 



s-1 



(1-*) 



Ff = 



f-s/2 



lb 



7i(ifo + C/i52 )(s-l) 



f1 
" 2 



_2 + s 



^oliKo 



s-^ — s 
2 



-(1-* 



,s + l 



1- 



(109) 

(110) 
(111) 
(112) 



^2 



1 



(1-*)^ 



(104) 



Since the frequency regime 1/4 < / < 1/4^1 is influenced 
by the solution of the part ^P^ < ^ < vl'g, for which we 
could not find an analytical expression, we negl ected the 
solutions in the overall result of the fiuence (Eq. (109l). 
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Nonetheless, the obtained result is notewortliy. As a matter 
of fact, it is remarkably simil ar to the case with s > 3 and 
ao ^ 1. The fluence of Eq. (109) is also represented by a 
broken power-law, with two different but constant spectral 
indices /3i = 1/2 and P2 = 3/2 for frequencies below unity, 
indicating that the non-linear cooling is replaced by the 
linear cooling for very late times. The spectrum is cut off 
at/ 



(7I, similar to every calculated case before. 



6.6. Numerical solution for a large injection parameter 

As we did for a <C 1, we want to compare our analytical 
result for ao 3> 1 with a simple numerical calculation. For 
the plots we dropped the constant Fa, which means that 
we normalized the fluence by NFqom — Fcom I Fa- 



a„ = 10 




12 3 

log(f) 

Fig. 22. NFcoM as a function of / for three cases of s 
{Full: s = 1.5, Dashed: s — 2.5, Dot-dashed: s — 3.5) and 
for ao = 10, and 52 = 10'^. The red lines indicate the an- 



alytical solution Eqs. (95 1 and (109 1, but with an offset of 
10-*. 



The results shown in Figs. [22] and [23] match each other 
rather well for / < 1 . The point / = Uq , where the power- 
law breaks from the linear to the non-linear cooling, is 
visible in the numerical and the analytical plots. However, 
in the numerical plot the spectral index for / < Oq is a 
little higher than 1/2. 

For / > 1 and s > 2 the numerical plots seem to follow 
an f^'^^^ dependence, and one should note especially that 
for s = 2.5 the plot is not as steep for / > 1 as for / < 1. 
This is expected, when we assume that the spectral index 
changes at / = 1 like A/32 — {s ~ 3)/2. 
For s < 2 the curve shows a similar behavior as in the 
case ap ^ 1, steepening gradually for / > 1, although 
the steepening is not as obvious as in the other case (cf. 
section 6.3 1. One should also note the small pile-up for 
s — 1.5 shortly before the cut-off, which occurs at f — g^ 
in all cases. 




At last we should note that the results of section |6| match 
those of SBM, at least for / < 1. This is reasonable, since 



J| -3 -2 -1 1 2 3 4 6 E 7 

log(f) 

Fig. 23. NFcoM as a function of / for three cases of s 
(Full: s — 1.5, Dashed: s = 2.5, Dot-dashed: s = 3.5) and 
for ao — 100, and 32 = 10'^. T he re d lines indicate the 
analytical solution Eqs. (95) and (109), but with an offset 
of 10-4. 

the power-law electron distribution quenches rather rapidly 
until it becomes a (5-function, which was used by SBM. 

7. Discussion and conclusions 

In this paper we have made further progress in the work 
on nonlinear cooling of highly relativistic electrons, as it 
could take place in jets of active galactic nuclei. We used 
the nonlinear approach of RS, who introduced a cooling 
scenario based upon the synchrotron self-Compton effect. 
We expanded that paper by using a power-law for the ini- 
tial electron distribution instead of a (5-function. 
We treated both the linear and the nonlinear cooling case, 
and calculated the respective electron number densities in 
section [2] which differed by a new time variable that en- 
tered the nonlinear solution. We showed that the nonlin- 
ear solution cools much faster for harder spectra than for 
softer ones. The linear solution also cooled rather rapidly. 
We discussed in section |2.3| that it depends critically on 
the choice of some parameters, which cooling acts quicker. 
This implies that the approach treating the cooling sce- 
narios separately is not necessarily fulfilled, because at one 
moment in time the linear cooling will always dominate 
over the nonlinear case. 

Therefore, we applied the electron distribution to a cool- 
ing scenario, which combines the linear and non- linear ap- 
proach in one term. This was also already calculated for 
a (5-distribution of electrons by SBM. We calculated the 
electron number density in section [3] and discussed that 
depending on the injection parameter ao the electron num- 
ber density showed characteristics of linear or non-linear 
cooling. However, even if the cooling is non- linear in the 
beginning, it becomes linear at later times, just as we ex- 
pected. 

We calculated the total synchrotron fluence for all cases and 
began with the linear and non-linear cases in sections [4[ and 
[5] respectively. In both cases the low-frequency regime of 
the fluence spectrum (/ < 1) does not depend on the in- 
jection spectral index s of the electrons. The fluence spec- 
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tral index is asYN = 1/2 for linear, and assc =^ 3/2 for 
nonlinear cooling, respectively. These constant spectral in- 
dices match those of the (5-function approach for the elec- 
tron number density used by RS. This is reasonable, since 
we showed in section 12.31 that the electron distribution be- 
comes a (5- function for late times, when the electrons have 
cooled significantly and the radiated synchrotron photon 
energy is low. 

For linear as well as for nonlinear cooling the spectra are 
cut off at / = 5|. In the region 1 < / < gf the spectral 
index is s/2 with the exception of the nonlinear case with 
small spectral indices, where the spectral index remains ba- 
sically unchanged from the low-energetic part. 
Section |6] was devoted to the calculation of the total syn- 
chrotron fluence of combined cooled electrons. The results 
are in some ways similar to the cases before. For ao <^ 1, 
which represents some kind of a linear solution, the flu- 
ence shows linear characteristics. Especially, below / = 1 
the spectral index is /3i = asYN — 1/2. For ao ^ 1, the 
fluence below / = 1 exhibits a power-law with a spectral 
index ^2 = o:ssc — 3/2, just as one would expect from 
a non-linear cooling. However, since the non-linear cooling 
becomes linear at later times, for frequencies below / — 
Oq" <C 1 the spectral index becomes /?i — osyn — 1/2, 
again. 

Our results prove that the spectral behavior of the total 
synchrotron fluences in the non-linear SSC and combined 
synchrotron-SSC cooling cases is practically independent 
from the functional form of the energy injection spectrum. 
The earlier predictions of SBM, based on a mono-energetic 
injection distribution functions of electrons, on the syn- 
chrotron fluence behavior therefore hold in this case. 
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Appendix A: Calculation of the non-linear electron 
number density 

Inserting the nonlinear cooling rate Eq. ^ and the injec- 
tion rate into Eq. ^ gives us the differential equation for 
the SSC cooled electron number density nssc (we drop the 
subscript in the following): 



dn{j,t) d 



dt 



dj 



Aoi' 



d77^n(7,t) 



nil,t) 



= qo-/-' H [7 - 7i] H [72 - 7] S (t) . (A.l) 

Multiplying the equation with j"^ /Aq we obtain with the 
definitions y = AqI and S 



7^n 



as; 

dy 



d-fS 



,dS_ 



Qol 



S [y) H [72 - 7] H [7 - 71] 



(A.2) 



We yield with ^ = 7 ^ 

dS_ 
dy 






as 



qoe-'Siy) H[C-6] H[a-C] 



If we define the implicit time variable T through 



dr 

dy 



de^> 



(A.4) 



the differential equation becomes 

^ + ^^- 10^-' S {y) H [e - 6] H [a - i] . (A.5) 

Formally multiplying this equation with dy / dT results in 
dS dS 



dT di 



qoe-'s{T) HK-e2] HKi-e] 



(A.6) 



This differential equation for the electron number density 
can be solved with the method of characteristics. Thus , we 
obtain 



d5 
dT 



X H Ki - T - Co] , 



(A.7) 



where ^0 = S, — T is a, constant of integration. Eq. (A.7 1 can 
be easily integrated with respect to T, which results in 

S - qo{^ - Ty-^ H [e - T - 6] H Ki - e + T] H [T] 

+Sh . (A.8) 

We now require that S'(^ — 0,T) — 0, which means 

Shi-T) = -qoi-ry-^ H [-T - 6] H Ki + T] H [T] (A.9) 

Collecting terms, we find S to be 

Sit T) = qoi^ - ry-' H K - T - 6] H [6 - e + T] 

x(H[T]- H[r-e]) .(A.IO) 

Since our fiare begins at T = 0, we are not interested in 
events that take place before that moment. Hence, we find 
the electron number density: 



n(7, T > 0) 



907^^(1-77^)'"' H 



X H 



72 



1 + 72T 



■7 



H 



7- 



1 
T 

71 



S 

■7 



I + 71T 



(A.ll) 



We see that we have two Heaviside functions defining upper 
limits for 7. It is an easy task to compare them, and to find 
out which one is lower than the other one. Having done 
so, we find the solution for the nonlinearly cooled electron 
number density to be 



"(7, 7^)= 907^^1 -7^)' ' H 

X H 



72 



I + 72T 

71 



I + 71T, 



(A. 12) 



(A.3) yielding Eq. ^. 
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Appendix B: Derivation of U 

The time variable T has been defined through Eq. (Isl): 
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U{T) 



dT 
dy 



d^TI 



(B.l) 



Inserting Eq. (A. 10 1, we gain 



oo 

u{T) = J d^r'qoi^ - ry-' H K - T - 6] 



X H [Ci + T - e] H [T] 



(B.2) 



As stated before, we are only interested in solutions for 
T > 0. Hence, we can neglect the third Heaviside function, 
which results in 



t/(T) = go J dee-2(e-r)^-^ 
A first substitution w — £, — T yields 

while a second substitution w ~ Tv gives 



(B.3) 



(B.4) 



6/T 



(7iT)-i 



qoT' 



dv 



(72T)-1 



{i + vy 



(B.5) 



where we re-substituted ^i — 7"^ in the last step. 
The purpose of the next two substitutions is to get rid of 
the time variable in the limits of the integral. In order to 
achieve this we first set a; = 71T and introduce 52 = 72/71, 
which yields 



3— s^s— 3 



U{x) = qolt "x' 



,s-2 



dv 



(g^x)- 



{l + vY 



l(^,,\s~2 



3-s [ A ^ ^^^^ 

Now, we use u = vx, resulting in 

1 



(B.6) 



U{x) = qo^i " J du ^^^^ 



-1 
92 



qolt" J du 



,s-2 



(x + uf 



(B.7) 



92 



This integral can be expressed in terms of the hypergeomet- 
ric function, but that would not yield an analytical form. 



Nonetheless, one can obtain an approximate solution in the 
regimes Q < x < g2 (small x), and x > 1 (large x). An 
analytical continuation serves as a solution for the interme- 
diate regime. For small x the integral can be written as 



Ui = (?o7i' 



,s-2 



du 



{x + uY 



92 



qolt' 



92 



9o7i 



,,3-s 



[1 



duu 



9t'] 



s-4 



s-3 

Similarly, we achieve for large x 

1 



(B., 



Us = qoii 



3-s 



du 



w 



5-2 



(x + uy 



qolt' 



92 



{S - 1).t2 



[1-52^-1 



(B.9) 



The requirement for the solution of the intermediate regime 
is that it must be continuous, meaning Ui{g2 ) = [^2(52" )i 
and C/2(l) = U-i{l). In order to accomplish such a behav- 
ior, we can first assume a proper solution U2 with some 
unspecified constants, and then try to fit it to the bound- 
ary conditions. A good ansatz is 



U2 



qoil " 



'^s-3 



h'x 



,92' 



(B.IO) 



Matching the solution with the boundary conditions, yields 
the values of the constants a' , b', and c': 



b' = 



a' = 1 
2 



s- 1 
s-3 



s-1 

Since we had only two equations for three parameters, we 

chose a' = 1. 

Thus, the obtained solution for the intermediate x-range is 



U. 



go7i ' 
s-3 



1 



2x 



s-3 



3 92 



s - 1 s-1 



(B.ll) 



Before we summarize the results, we need to say a few words 
about the spectral index s. We already stated that it must 
be greater than 1. But according to our results above, we 
also find that s ^ 3. Thus, we have two different cases to 
consider: s > 3, and 1 < s < 3. 
Collecting terms, we find for s > 3 




s-3 92 

s— 1 2;^ 



{s-l)x- 



[1-52^-1 



C7(X): 

0<x< g^' 
92' <x<l 
x> I 



(B.12) 
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and for 1 < s < 3 



s-> 3. 

For the third case we find 




U{x)-- 
, x>l 



yix>l) 



dx 



(B.13) 



7iC/0(,^(l-5r^) 
s-1 



37iC/o(s-3)(l-gr') 



C3 



(C.3) 



s-3 



and Ui — 



3-s 



Appendix C: The non-linear time variable x 

Since U{x) — '^^ i "^^ can separate the variables obtaining _ 

dy = — in^- This can be integrated quite easily except for 
both intermediate cases. However, we can find approxima- 
tive solutions by using the same approximations of these 2/2 — 
cases we used already during the calculation of the syn- 
chrotron spectra. 
Beginning with the case s > 3 we yield 



These equations can be inverted simply, yielding x{y). As 
for U{x), we require x to be continuous at the points yi = 
y(x — §2^) and y2 = y{x — 1). Matching the solutions at 
these points we find the values 



.92 



y{0<x< 52-1) = 



dx 



7iC/o(i - gtl 



1 



7iC/o(l"52~') 



-X + Ci 



(C.l) 



We require y{x = 0) = 0, which means ci = 0. The inter- 
mediate range becomes 



2/(ff2"' < a^ < 1) = 



dx 

7iC/o ( 1 - '"''"' 



.s-l 



1 



liUo 



dx — 



liUo 



s-3 92 
s— 1 x^ 



X + C2 



(C.2) 



5" 



■ s =3.1 

3 =4 
- s = 5 




-4 -3.S -3 -2.6 -2 -1.6 -1 -0.6 

log{x) 

Fig. C.l. The denominator U{x) of the intermediate time 
regime as a function of x for three cases of s. 



C2 



C3 



1 

liUo 

1 


92 

r i-.92'-^J 

91- 


71^0 ] 

1 
71^0 


~ 9t' 

[ l-9t' 



s- 1 



3{3-s){l-gl-n 
For the case 1 < s < 3 we find similarly 

y{0<x< 32-1) ^ 



(C.4) 
(C.5) 
(C.6) 

(C.7) 



7iC/i(.g2'-^-l) 
1 

X + Ci . 



7iC/i(52''^-l) 



(C.8) 



We require y{x = 0) = again, which means Ci = 0. The 
intermediate range becomes 



2/(52^' < X < 1) 
1 



dx 



^, Tj ( 2x^zl 



s-3 9-2 

s-l x'^ 



llUo 



dx 



„4-s 



27iC/i(4-s) 



C5 



(C.9) 



5" 













X^^ 










10= 


■^X. 


* s = 2 
s = 5/2 


- 


10^ 




in^ 






~'~^-- **•» ^^\- 




~^^^^--==^.^^l*;-*^^x. 














111 


"^ --. 








1 1 1 1 1 1 1 



-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

log(x) 



We approximated the denominator with the leading term. Fig. C.2. The denominator U{x) (black) and its leading 
The validity of this approximation for most parts of the term (red) of the intermediate time regime as a function of 
x-range can be seen in Fig. |C.1[ It becomes less valid for x for three cases of s. 
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We approximated the denominator witli tlie leading term. D.l. Large spectral index 

The validity of this approximation for most parts of the x- „. ., , . i- r^ 

u • !?• iTT-oi rjiT, 1 i r Similarly to the steps m appendix C we calculate the de- 

ranee can be seen m i-iff. L. 2 i he largest errors occur tor , -/ x ,tt , • •, , ' — ' „ t^ ~ -i 

small values of x. 



The last case yields 
y(x > 1) = 



pendence x{y). We begin with the case s > 3. For x < §2 
we find 



dx 



71 c/i (3^(1 -.9r) 

s - 1 

-C6 . 



y{x < 32 ^) = 



1 

7i 



d.T 



37it/i(3-s)(l-.9r') 



(CIO) 



Ko + C/o(l - gt') 

X 

■ Ci . 



MKo + Uoil-gtn) 



(D.5) 



As before, a simple inversion leads to x{y), while the re- 
quirement that the solution should be continuous at the 



As b efore, we set Ci = 0, since y{x = 0) = 0. Inverting Eq. 
dRSl) yields 



points 2/3 = y{x = g^ ^) and y^ = y{x = 1) gives the values x{y < yi) = 7i(ii:o + Uo{l - gl^''))y 



J/3 



Vi 



C5 



ce 



1 92' 



7if/i gr^- - 1 


1 


■ s-1 


71C/1 


[2(4 -s) 


1 


r .92-^ 


71C/1 


U- 1 


1 


■ s-l 



(1-52 ) + ^3 



.92 



92 



3-s 



2(4- 



■s).92~ 



71C/1 



2(4 -s) 
s-1 



(i-5r*) + 



52-^ 



52^-^ 



3(3-,s)(l-.9r'^) 



(C.ll) 
(C.12) 
(C.13) 

(C.14) 



Obviously, yi is found from the condition x[y 
yielding 

y.^g^'hiiKo + Uoil-gtl))-' ■ 
For g2 < a: < 1 we find 

1 r dS: 

71 J Ko + Uo 

X 

C2 , 



yi) 



(D.6) 
-92' 

(D.7) 



2/(52^' < i < 1) = 



Appendix D: The implicit time variable x of the 
combined cooling 

The differential equation (pi may look a little bit more com- 
plicated with the combinea cooling term ( 13 ). However, the 
solution can be obtained with the methods outlined in ap- 



or the other way around 

x{yi < 2/ < 2/2) = 7i(-f^o + Uo){y - £2) 



(D., 



(D.9) 



Since x is supposed to be continuous, we find the constant 
C2 by matching the solutions for y — yi resulting in 



C2 



U09I-' 



pendix 



A yielding solution (14 1. The important difference 



is the dehnition of the implicit time variable, which has to 
be chosen as 



jiiKo + Uoil-gtlKKo + Uo) ' 
We also obtain 2/2 from the condition x{y — y2) — 1 

1 /Xo + C/o(l 



U 



dr 
dy 



2/2 



Ko+ djj^ncoMil.r) 



(D.l) 



.92^-^ 



.92^-^^ 



Using the definition of U from appendix |B] this can be 
written as 

U ^Ko + U . (D.2) 

Thus, we can use the previous results to obtain for s > 3 

U{i,s > 3) = 



li{Ko + Uo)\ Ko + Uo{l-gtn 
Defining Bq = f/ofiif (1 - .92 '*) we yield for i > 1 
/ - .-. 1 /" dx 

y{x > 1) 



(D.IO) 



(D.ll) 



Ko + Ua [I 
Ko + Uo 

Kn + U, 



9l-'] 



and for 1 < s < 3 



[1-.92^~1 



< X < g^' 

52"' < i < 1 
x> 1 



liKo 



■ arctan 



71 7 Ko + So£-2 


Wdi ^' 


£2 


7ii^o y ""^ t + 







(D.12) 




U{x, 1 < s < 3) = 

,0<i<g2^ 

, 92^ < i < 1 
, X > 1 



(D.3j T^lie problem arising is that we cannot find an inverted ex- 
pression for X. However, we can obtain approximative re- 
sults for small and large arguments of the arctan-function. 
In order to achieve these approximations we define the in- 
jection parameter 



(D.4) "0 




907f"'(l 



92 ) 



Kois - 1) 



where we defined x = 71 r, and used the leading term ap- 
proximation discussed in appendix [O for the intermediate 
regimes. 



for which Eq. (D.12| becomes 
1 



y{x > 1) 



71-f'^o 



ao arctan 



+ ^3 



(D.13) 



(D.14) 
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For ao <^ I the argument of the arctan-function is always ao ^ 1 the solution is purely linear, while for ao 3> 1 it 

(much) larger than unity, since x > 1. We, therefore, ap- is non-linear and becomes linear at later times, just as we 

proximate arctan(a;/ao) ~ 7''/2, set 63 = £3 — 7r/2, and expected it. 

obtain Before we proceed with the case 1 < s < 3 we list the re- 

1 suits of this section once more in a compact form. We also 

y{x > l,ao ^ 1) = — 71- [x + £3] . (D.15) approximate the results for 32 ^ 1, which, as one will see, 



71^0 
This can be easily inverted, yielding the linear solution 

x{y>y2,aa'^^)=liKa{y~C4)- (D.16) 



Matching this solution with Eq. (|D.9|) yields 
1 



C4 = 



ji{Ko + Uo) 



Ko + Uoil + gt'~92) 



Ko + Uo{l-gtn 



1 



71-^0 



(D.17) 



simplifies a lot. 

x{y < yi) « 7i{Ko + Uo)y 

x{yi <y<y2) = ji{Ko + Uo){y - C2) 

^(y > y2,ao < 1) = liKo{y-C4) 

xim < 2/ < y3,ao > 1) = [3ao7i^o(2/- C5)] 

x{y > y3,ao > 1) = jiKoiy-ce) , 
with 



For ao ^ 1 we have to consider two cases. If 1 < x < ap, we 
see that x/ao < 1. Thus, we can approximate the arctan- 
function to third order as arctan(a;/ao) ~ x/aQ — x'^/Sag, 
resulting in 



1 



JiKo 



X — ao 



y(l < X < ao,ao > 1) 
X x"^ 
"0 3ag 
1 

2 +C3 



C3 



71^0 



3a 







(D.18) 



Inverting yields 

xim < y < 2/3, ao > 1) = [^alliKoiy - C5) 
with 

"0 , 1 



2/1 R 
h - 

^3 R 
and 

C2 ft 
C4 R 



92 

7i(A'o + C/o) 
1 

1 



al-l 



1/3 



(D.19) c~5 



C6 



ji{Ko + Uo) 


ialjiKo 


UogI-' 


5. n 


7i(Xo + f/o)2 




Uo 




-fiKo{Ko + Uq) 


(3a2 - l)Ko 


-Uo 


3al-fiKo{Ko + Uo) 


1 


2al + l 



2/3 



371 i^o 7i(^"^'o + f^o) 



-/i{Ko + Uo) ial-iiKo 



(D.25) 
(D.26) 
(D.27) 

(D.28) 
(D.29) 

(D.30) 
(D.31) 
(D.32) 

(D.33) 
(D.34) 
(D.35) 
(D.36) 



^0 + ^:^0(1+272-52"') 



1 



and 

C5 = 



Ko + Uo{l-9l-') ) ial-tiKo' 
1 /i^o + C/o(l + 2/2"'-272"') 



(D.20) 



71(^0 + C/o) 



i^o + C/o(l-2/2 ) 



1 



iaf^-iiKo 



(D.21) 



Since in this approximation C2 — >■ 0, x{y < yi) = x{yi < 
y ^ 2/2), and, thus, one can neglect yi. 

D.2. Small spectral index 

We will now derive the explicit form of the implicit time 

variable x for 1 < s < 3. 

The first regime is x < 5^ , yielding 

,- -In 1 /" dx 

y[x < 32 ) = 



In the case x > ao, 'we can approximate again 
arctan(x/Q;o) ~ 7i'/2, yielding with £3 = £3 — tt/2 



11 J Ko + Ui{gl-^-l) 

X 



2/(x > ao,ao > 1) 



1 



li{Ko + UM'' 



(D.37) 



[^ + 23] , 



7i-f^o 
or inverted 

S:{y> y3,ao > 1) ^ -yiKoiy ~ ce) , 
with 

1 fKo + Uo{l + 9l-'-gtn 



2-1)) 

(D.zzj Since 2/(2: = 0) = 0, obviously di = 0, and the inversion 
becomes 



(D.23) 



xiy <2/i) = 71(^0 + C^i(ffi 
with 



3-s 



1))?/ 



C6 



71(^0 + C/o) 



ifo + C/o(l-.92'"') 

1 2ao 



yi=7i(i^o + [/i (52^^-1)). 

The next time step is 5^ < x < 1, resulting in 



(D.38) 
(D.39) 



Sa^jiKo 3-fiKo 



/ -1 ~ 1 \ 1 /" dx 

(D.24) 2/fe <x<l) = -y^^-^^ 



s-3 



What one can see here is that the injection parameter con- 
trols significantly the cooling behavior of the electrons. For 



271 C/i 



dx 



A + x'-^ 



(D.40) 
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where we defined /3 = (3 — s)(l — g\^'')/2. For ao ^ 1 we where 
see that fi/a^ ^ 1 (as long as s is not too close to 1 or 3) 
and with a; < 1 we approximate the integral as 



y{g2 ^ < 5 < l,ao < 1) 



^ s-1 
"" 271(71 J /3 
(s - 1)^2 



(s-l)a| _i. 



.92 



Odx 

y + ^2 • 



7i(i^o + C/i(.92'"^-l)) 
1 



7i-f''o 



(D.52) 



27iC^i/5 
The inversion is easily performed, yielding 

271 C/i^ , 
x{yi < y < y2:"o< 1 = 7 t^^ y-a2 

(s- l)a^ 
where we obtain by matching the solutions 

5/2 = o, rr Q (l^g2 ) 



(D.41) 



(D.42) 



For ao 3> 1 we use for 1 < a; < ao the approximation 
arctan(i/ao) ~ xjao — a;'^/3Q:Q, yielding 

1 



2/(1 < a; < ao,Q:o > 1) 



Hence, 



71^0 



^.^'^ 



y2 



and 



d2 = 



92 



271 C/i/? 



.92 



7i(ifo + c/i(.g2'"' 

(s - l)ai 



1)) 



,(D.43) ^(2^3 <2/ < y4>"o > 1) =^ [37iii'oao(y-rf6)] 
with 

>3-i {s-m-gr') 



1/3 



(D.53) 



(D.54) 



(D.44) 



71 V^o + (7i(.g3-^-l) 2t/i/3 

For ao ^ 1 we see that /^/ofg <C 1. As a rough approxima- 
tion this is also much lower than 172^ , and, therefore, we 
achieve the integral 



Va 



and 



■^0-1 

SjiKnal ' 27iC/i(4-s) 
-1 



.92 



2/(32 ^ < i < l,ao > 1) 



s-1 



2jiUi 



^s-3 1 - 



x" ■^ dx 



71(^0 + (7i(.9r'-l)) 



(.-l)(l-,gr') 



(D.55) 



.92-^ 



s-1 



27iC/i(4-s) 7i(ifo + C/i(.92'"^-l)) 



;^4-s 



27iC/i(4-s) 
The inverted equation is 

27i;7i(4-s) 



xivl < 2/< 2/3>ao > 1) 
with 

(i-5r') 



s- 1 



(2/-rf3) 



(D.45) 



(D,46) 



1 



3^1 Koa^ 



(D.56) 



The last case is for x > aQ, where we can use the "linear" 
approximation arctan(a;/Q;o) ~ 7''/2, again. With d"' = d^ — 
7r/2 we achieve 



y{x > ao,ao > 1) 



1 



s-1 



71^0 



li + d'l 



and 

d.^ 



27iC/i(4-s) 
92' 



-tiiKo + Urigt"-!)) 



Jh 

,3-s 



or inverted 
(D.47) x(2/ > yl,ao > 1) - 71^0(2/ - ^7) , 



(D.57) 



(D.58) 



where we defined 



(.'^-1)2/2 



s-4 



7i(i^o + Uiigl-' - 1)) 271(71 (4 - s) 



(D.48) ^7 - 



92 



Similarly to the case i > 1 for large spectral indices, we 
obtain here for the integral in that time regime 



(s-l)(l-gr") 

27it/i(4-s) 71 (/fo + C/i (.92^^-1)) 

2al + l 



S-yiKoa^ 



(D.59) 



y{£ > 1) 



1 



x'^ dx 



jiKo J al + x"^ 



1\Kq 



X — a^ arctan 



an 



As we did for the case of large spectral indices, we sum up 
our results in a short list, and perform the approximation 
for 52 > 1. 



(D.49) 



i{y<yl) ~ MKo + Uig^ ')y 
271C/1/3 



We will continue with the same approximations as before, ^(wi* < w < yo ao <^ 1) = ~ '"^" "^^ (y — ^2) 
-..-., s „,.,,„ - - , (s-l)ag 



yielding for ao ^ 1 arctan(a;/ao) ~ 7r/2, and with d'^ — 
di — 7r/2 the result 



1 



-(x + 4) 



y{x > 1,q;o < 1) „ , 

71^0 

The inversion is obviously 

x{y > y*2,aa < 1) = jiK^iy - d^) , 



x{y >y*2,aQ<^l) = 'fiKoiy - dr,) 
27i;7i(4-s) 



(D.50) x(yi*<2;<y3*,ao»l) 



s- 1 



(2/ - ^3) 



^(2/3 < 2/ < y4>"o > 1) = [37i^oao(y-c?6)] 
(D.51) x{y>yl,ao:$>l)^'fiKo{y-dj), 



(D.60) 
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(D.62) 
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(D.64) 
(D.65) 



24 

with 

~* -^ , 
Vi ~ 

J/3 ~ 

and 

d2 « 
(^5 ~ 

dr - 



92 



liiK^ + Uagt 
{s-l)al 



3-s\ 



.92 



271C/1/3 7i(i^o + C/og2 ) 



s- 1 
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D.66) 
D.67) 
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